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Bohmian trajectories from a Dynamical Systems point of 



Abstract 

Vortices are known to play a key role in the dynamics of the quantum trajectories de- 
fined within the framework of the de Broglie-Bohm formalism of quantum mechanics. It 
has been rigourously proved that the motion of a vortex in the associated velocity field can 
induce chaos in these trajectories, and numerical studies have explored the rich variety of 
behaviors that due to their influence can be observed. In this paper, we go one step fur- 
ther and show how the theory of dynamical systems can be used to construct a general and 
systematic classification of such dynamical behaviors. This should contribute to establish 
some firm grounds on which the studies on the intrinsic stochasticity of Bohm's quantum 
trajectories can be based. An application to the two dimensional isotropic harmonic oscil- 
lator is presented as an illustration. 
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1 Introduction 

Some interpretational difficulties [1] with the standard version [0 led David Bohm to develop 
in the 1950's (3]| an alternative formulation of quantum mechanics. Despite initial criticisms, 
this theory has recently received much attention 0] |5]|, having experimented in the past few 
years an important revitalization, supported by a new computationally oriented point of view. 
In this way, many interesting practical applications, including the analysis of the tunnelling 
mechanism J6l|71[8]|, scattering processes fl9l [TOl HTJ , or the classical-quantum correspondence 
lfT2l[T3ll , just to name a few, have been revisited using this novel point of view. Also, the chaotic 
properties of these trajectories lfT4l [T51 [T6l [P71 [T8l , or more fundamental issues, such as the 
extension to quantum field theory |fT9l , or the dynamical origin of Born's probability rule ll20l 
(one of the most fundamental cornerstones of the quantum theory) have been addressed within 
this framework . 

Most interesting in Bohmian mechanics is the fact that this theory is based on quantum 
trajectories, "piloted" by the de Broglie's wave which creates a (quantum) potential term addi- 
tional to the physical one derived from the actual forces existing in the system [[3J . This term 
brings into the theory interpretative capabilities in terms of intuitive concepts and ideas, which 
are naturally deduced due to fact that quantum trajectories provide causal connections between 
physical events well defined in configuration and time. Once this ideas have been established 
as the basis of many numerical studies, it becomes, in our opinion, of great importance to pro- 
vide firm dynamical grounds that can support the arguments based on quantum trajectories. For 
example, it has been recently discussed that the chaotic properties of quantum trajectories are 
critical for a deep understanding of Born's probability quantum postulate, considering it as an 
emergent property ll20ll . Unfortunately very little progress, i.e. rigorous formally proved math- 
ematical results, has taken place along this line due to the lack of a solid theory that can foster 
this possibility. Moreover, there are cases in the literature clearly demonstrating the dangers of 
not proceeding in this way. One example can be found in Ref. GT1 . where a chaotic character 
was ascribed to quantum trajectories for the quartic potential, supporting the argument solely 
on the fact that numerically computed neighboring pairs separate exponentially. This analy- 
sis was clearly done in a way in which the relative importance of the quantum effects could 
not be gauged. Something even worse happened with the results reported in Il22ll . that were 
subsequently proved to be wrong in a careful analysis of the trajectories ||23l . 

Recently, some of the authors have made in Refs. lfT5l [T6l [T71 what we consider a relevant 
advance along the line proposed in this paper, by considering the relationship between the even- 
tual chaotic nature of quantum trajectories and the vortices existing in the associated velocity 
field which is given by the quantum potential, a possibility that had been pointed out previously 
by Frisk [TT4l . Vortices has always attracted the interest of scientists from many different fields. 
They are associated to singularities at which certain mathematical properties become infinity or 
abruptly change, and play a central role to explain many interesting phenomena both in clas- 
sical and quantum physics [|25l . In these papers it was shown that quantum trajectories are, in 
general, intrinsically chaotic, being the motion of the velocity field vortices a sufficient mecha- 
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nism to induce this complexity |fT5l . In this way, the presence of a single moving vortex, in an 
otherwise classically integrable system, is enough to make quantum trajectories chaotic. When 
two or few vortices exist, the interaction among them may end up in the annihilation or creation 
of them in pairs with opposite vorticities. These phenomenon makes that the size of the regular 
regions in phase space grows as vortices disappear IfTTTl . Finally, it has been shown that when 
a great number of vortices are present the previous conclusions also hold, and they statistically 
combine in such a way that they can be related with a suitably defined Lyapunov exponent, as a 
global numerical indicator of chaos in the quantum trajectories lfT6ll . Summarizing, this makes 
of chaos the general dynamical scenario for quantum trajectories, and this is due to the existence 
and motion of the vortices of the associated velocity field. 

In this paper, we extend and rigorously justify the numerical results in |fT6l[T7l[T5l concern- 
ing the behavior of quantum trajectories and its structure by presenting the general analysis of 
a particular problem of general interest, namely a two-dimensional harmonic oscillator, where 
chaos does not arise from classical reasons. In this way, we provide a systematic classification 
of all possible dynamical behaviors of the existing quantum trajectories, based on the appli- 
cation of dynamical systems theory [|24l . This classification provides a complete "road-map" 
which makes possible a deep understanding, put on firm grounds, of the dynamical structure for 
the problem being addressed. 



2 Bohmian mechanics and quantum trajectories 

The Bohmian mechanics formalism of quantum trajectories starts from the suggestion made by 
Madelung of writing the wave function in polar form 

il>(r,t) = i?(r,t)e i5 ^, 

where R 2 = ipijj and S = (Imjj — ln-0)/(2i) are two real functions of position and time. 
For convenience, we set h = 1 throughout the paper, and consider a particle of unit mass. 
Substitution of this expression into the time-dependent Schrodinger equation allows to recast 
the quantum theory into a "hydrodynamical" form Q, which is governed by, 

^-v.(*vs), (1) 
as jvsy _ iv 2 r 

dt ~ 2 V 2 R ' { ) 

which are the continuity and the "quantum" Hamilton- Jacobi equations, respectively. The qual- 
ifying term in the last expression is customarily included since this equation contains an extra 
non-local contribution (determined by the quantum state), Q = ^V 2 R/R, called the "quantum" 
potential. Together with V, this additional term determines the total forces acting on the system, 
and it is responsible for the so-called quantum effects in the dynamics of the system. 
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Figure 1 : Illustration of the dynamical consequences of a time-reversible symmetry. 



Similarly to what happens in the standard Hamilton- Jacobi theory, Eqs. © and © allow 
to define, for spinless particles, quantum trajectories by integration of the differential equations 
system: f = — W(r) — VQ(r). Alternatively, one can consider the velocity vector field 



X^ = VS = - Y Y ]L J Y . (3) 



1 ipWip — i/jWi/j 

2 W 

Notice that, in general, this Bohmian vector field is not Hamiltonian, but it may nevertheless 
have some interesting properties. In particular, for the example considered in this paper it 
will be shown that it is time-reversible, this symmetry allowing the study of its dynamics in a 
systematic way. 

Let us recall that a system, r = X(r,t), is time-reversible if there exists an involution, 
r = 0(s), that is a change of variables satisfying O 2 = Id and 9 ^ Id, such that the new 
system results in s = DO~ 1 (s)X(B(s),t) = —X(s,t). One of the dynamical consequences 
of reversibility is that if r(t) is a solution, then so it is 0(r(— t)). This fact introduces sym- 
metries in the system giving rise to relevant dynamical constraints. For example, let us as- 
sume that 6(x, y) = (x, —y) is a time-reversible symmetry (see Fig. [[])• Then any solution 
r{t) = (x(t),y(t)) defines another solution given by (x(— t), — y(— t)). Let us remark that 
this fact constraints the system dynamics since if, for example, r{t) crosses the symmetry axis 
(y = is invariant under 6) then the two solutions must coincide. 

We conclude this section by stressing that time-reversible systems generated a lot of interest 
during the 80's due to the fact that they exhibit most of the properties of Hamiltonian systems 
(see l|26l |27l |28l ). In particular, this type of systems can have quasi-periodic tori which are 
invariant under both the flow and the involution 0. That is, KAM theory fully applies in this 
context. Furthermore, some interesting results concerning the splitting of separatrices have been 
developed successfully for time-reversible systems ||29l , providing powerful tools for the study 
of homoclinic and heteroclinic chaos. 
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3 Model and canonical form 

The system that we choose to study is the two dimensional isotropic harmonic oscillator. With- 
out loss of generality, the corresponding Hamiltonian operator for r = (x, y) can be written in 
the form 

In this paper, we consider the particular combination of eigenstates: <^ ,o = 1/ with energy 
1, and (fifi = 2x/V2n, <^o,i = 2y/y/2n with energy 2. It can be immediately checked that the 
time evolution of the resulting wave function is given by 

, fAe- [t 2xBe~ 2it 2yCe~ 2it \ 

jP = + -=- + y e ^ x +y \ (4) 



whereQ A = A + iD, B = B + iE and C = F + iC, subject to the usual normalization condition 
\A\ 2 + \B\ 2 + |C| 2 = 1. In addition, we further assume the condition BC ^ EF in order 
to ensure the existence of a unique vortex in the velocity field at any time. Accordingly, the 
quantum trajectories associated to ® are solutions of the system of differential equations: 

-2(BC - EF)y - ^2(BD - AE) cost - s/2(AB + DE) sint 

x = T77 7^ ' (5) 

V{x,y,t) 

2{BC - EF)x + V2(AC - DF) cos t - V2(DC + AF) sin * 

, (6) 



V(x,y,t) 



where 



V{x,y,t) = 2{B 2 + E 2 )x 2 + 2{C 2 + F 2 )y 2 + A{BF + EC)xy + D 2 + A 2 
+2V2((AB + DE) cost + (AE - DB) sint)x 
+2v / 2((£>C + AF) cost + (AC - DF) sinf)y. 

To integrate this equation a 7/8-th order Runge-Kutta-Fehlberg method has been used. More- 
over, since the vector field is periodic, the dynamics can be well monitored by using strobo- 
scopic sections. In particular, we plot the solution (x(t),y(t)) at times t = 27m for n = 
1, 2, . . . , 10 4 and for several initial conditions. 

In Fig. [2] we show the results of two such stroboscopic sections. As can be seen the left plot 
corresponds to completely integrable motions, whereas in the right one sizeable chaotic zones 
coexisting with stability islands, this strongly suggesting the applicability of the KAM scenario. 
However, our vector field is neither Hamiltonian nor time-reversible, and then the KAM theory 
does not directly apply to this case. However, we will show how a suitable change of variables 



'The choice of notation for the real and imaginary parts of C may look arbitrary at this point, but it makes 
simpler the notation for the canonical form introduced in the next section. 
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Figure 2: Stroboscopic 27r-periodic sections for the 
different values of the normalized constants. 
Left plot: A = 0.37, D = -0.02, B = C = 0.44 and E 
Right plot: A = 0.4, D = -0.018, B = 0.37 and C = 1 
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quantum trajectories generated by Eqs. (0 and © for 

= —F. 
5 = 0.49. 



can be performed that unveils a time-reversible symmetry existing in our vector field. For this 
purpose, we first recall that the structure of gradient vector fields is preserved under orthogonal 
transformations. In this way, if we consider the transformation r = Ms, with M T = M _1 , 
applied to r = X7S(r, t), we have that s = VS(s, t), being S(s, t) = S(Ms, t). In other words, 
any orthogonal transformation can be performed on the wave function instead of on the vector 
field. 

Lemma 3.1. If Eq. © satisfies the non-degeneracy condition BC ^ EF, then there exist an 
orthogonal transformation and a time shift, such that the wave function takes the form 

, / ie-« 2xBe~ 2it 2yiCe~ 2[t \ .u^,^ 

V V 71 " V27T V27T / 

where A, B ,C G R, B > 0, C + satisfy A 2 + B 2 + C 2 = 1. 

We will refer to the wave function © as the canonical form of ©, and the rest of the paper 
is devoted to the study of this case. For this reason, the hat in the coefficients will be omitted, 
since it is understood that D = E = F = 0. In Table CD we give the actual values the canonical 
coefficients after the transformation corresponding to the results in Fig. [2l 

Proof. For convenience, we consider the complexified phase space z = x + iy, so that the wave 
function © results in 

if>(z,z,t) = (Ae^ + Be^z + Ce-^zAe-^, 
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Left plot 


Right plot 


A 


0.370540146272978 


0.400404795176082 


D 


0.656772411113622 


0.705788460189184 


C 


0.656772411113622 


0.584413081188110 



Table 1 : Wave function coefficients in the canonical model corresponding to the results of Fig. [2] Hats have 
been omitted as discussed in the text. 



where ^A = A + iD, \f^B = B + C + i(E — F), and yfh£ = B-C + i{E + F). Then, 
it is easy to check that the vortex, i.e. the set of points where the wave function vanishes, has 
the following position with respect to time 



z v (t) = 



-\A\\B 


e i(<- 


b+a) , 


lillcV^-^) 




B 


2 _ 


c 


2 



where A = \A\e m , B = \B\e lb and C — |C|e lc . Notice that the vortex is well defined thanks to the 
non-degeneracy assumption, and its trajector}]! follows an ellipse. This ellipse does not appear 
in the usual canonical form, but this can be made so by performing the rotation: z i— > ze 1 ^ and 
the time shift: t \— > t + A. In this way 



z v (t) = 



-\A\\B 


gi(t— b+a—. 




X) + 


|„4||C|e^'- c+a +^+ A ) 




B 


2 _ 


c 


2 



where it is clear that by choosing 2/i = c — b and 2A = c + b — 2a, the desired result is obtained. 
Then, the corresponding wave function in these new coordinates is 

^ = (\A\e- H + \B\e~ 2it z + ICIe-^e-^V^"^), 

that can be further simplified since the factor e^ 2a -^) plays no role in the Bohmian equations 
for the quantum trajectories. Finally, by recovering the coefficients in cartesian coordinates, one 
obtains A = y/n\A\, B = ^\{\B\ + |C|), C = (\B\ - \C\) and D = E = F = 0, which 
renders Eq. ©. □ 



4 Study of the canonical form 

Throughout the rest of the paper we consider the wave function © with A, C ^ 0, B > and 
D = E = F = 0. Let us remark that by changing the time t i— ► —t, if necessary, we can further 



2 We use here the term trajectory to refer to the evolution of the vortex, despite the fact that it is not a solution 
of the ODE. 
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restrict the study to the case C > 0. The corresponding quantum trajectories are then obtained 
from the vector field 



where V(x,y,t) = 2B 2 x 2 + 2C 2 y 2 + 2y/2ABxcost+2V2ACysmt+A 2 . In these coordinates, 
the only vortex of the system follows the trajectory given by 



which corresponds to an ellipse of semi-axes a = A/{y/2B) and b = A/{y/2C), respectively. 

In Fig. [3] we show some stroboscopic sections corresponding to this (canonical) velocity 
field for different values of the parameters a and b. As can be seen, a wide variety of dynamical 
behaviors, characteristics of a system with mixed dynamics, is found. In the left-top panel, 
which corresponds to the case in which a = b (vortex moving in a circle), we have sections 
corresponding to a totally integrable case. As we move from left to right and top to bottom 
some of these tori are broken, and these areas of stochasticity coexist with others in which the 
motion is regular, this including different chains of islands. Moreover, the size of the chaotic 
regions grows as the value of b separates from that of a. 

This variety of results can be well understood and rationalize by using some standard tech- 
niques of the field of dynamical systems, in the following way. Although the vector field © is 
not Hamiltonian, it is time-reversible with respect to the involution Q(x, y) = (x, —y). This re- 
sult is very important for the purpose of the present paper, since it implies that the KAM theory 
applies to our system if we are able to write down our vector field in the form X$ = X + eX 1 , 
£ « 1, being the dynamics corresponding to X integrable and Xi time-reversible. More 
specifically, let us assume that X does not depend on t and Xi be 27r-periodic with respect to 
t. Moreover, let us assume that for X there exists a family of periodic orbits whose frequency 
varies along the family (non-degeneracy condition). Then, our result guarantees that when the 
effect of the perturbation eXi is considered, most of the previous periodic orbits give rise to 
invariant tori of frequencies (1, uj), where uj is the frequency of the unperturbed periodic orbit. 
Of course, the persistence of these objects is conditioned to the fact that the vector (1, uj) satis- 
fies certain arithmetic conditions (see Il27ll28l for details). Since these arithmetic conditions are 
fulfilled for a big (in the sense of the Lebesgue measure) set of the initial orbits, the important 
hypothesis that we have to check in order to ascertain the applicability of the KAM theory is 
the non-degeneracy of the frequency map. 

In our problem, two such integrable cases exists. First, if A = the vortex is still at 
the origin and the time periodic part in the vector field disappears. As a consequence, all the 
quantum orbits of the system appear as ellipses centered at the origin in the xy-plane. It will 
be shown in the next section that the corresponding frequency varies monotonically along the 
orbits. This case has not been explicitly included in Fig. [3] due to its simplicity. Second, and 




(8) 
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Figure 3 : Stroboscopic 27r-periodic sections corresponding to the quantum trajectories generated by the canonical 
velocity field defined by Eq. ® for a = 0.4 and, b = 0.4, 0.44, 0.48 and 0.68 from left-top to right-bottom, 
respectively. 

as will be analyzed in Sect. [6l if B = C, or equivalently a = b, that is the vortex moves in 
a circle, the vector field is also integrable for any value of A. The corresponding stroboscopic 
sections are shown in the top-left panel of Fig. [3]). Here, the structure of the phase space changes 
noticeably, since two new periodic orbits, one stable and the other unstable, appear. Moreover, 
the obtained integrable vector field depends on t. We will show that this time dependence 
can be eliminated by means of a suitable change of coordinates, showing that our problem 
remains in the context described in the previous paragraph. The rest of the panels in Fig. [3] 
can be understood as the evolution of this structure as the perturbation, here represented by the 
difference between B and C, as dictated by the KAM theorem. 

To conclude the paper, let us now discuss in detail the two integrable cases in the next two 
sections. 
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5 The integrable autonomous case 

For A = 0, B ^ and C 7^ it is easily seen that the vector field © is integrable. Actually, 
the orbits of the quantum trajectories in the icy-plane are ellipses around the origin (position 
at which the vortex is fixed). Also, the frequency of the corresponding trajectories approaches 
infinity as they get closer to the vortex position. Let us now compute the frequency u of these 
solutions. First, we introduce a new time variable r, satisfying dt/dr = B 2 x 2 + C 2 y 2 , and then 
solve the resulting system, thus obtaining 

x(t) = acos{BCr + (3), y(r) = asm{BCr + (3), (9) 

where a is the distance from the vortex. Next, we recover the original time, t, by solving the 
differential equation defining the previous change of variables 

dt r > B 2 + C 2 9 B 2 — C 2 „ „ 

— = a 2 h a 2 cos(2BCt + 2/?), 

GST 2 2 

whose solution is given by 

2 B 2 - C 2 



a 2 (B 2 + C 2 ) 2(B 2 + C 2 )BC 



sin(25CV + 2/3). 



Notice that this equation is invertible since \2BC5\ < 1, and then r = jt + f(2BC n ft), f being 
a 2ix -periodic function. Finally, introducing this expression into ©, one can conclude that the 
solution has a frequency given by 

2BC 

U = BCl= a 2 {B 2 + C 2 ) 

that varies monotonically with respect to the distance to the vortex. Then, for A 1, the 
existence of invariant tori around the vortex is guaranteed. 



6 The integrable non-autonomous case 

Let us consider now the case of non- vanishing values of B = C for any A ^ 0. In this case, we 
have a = b 7^ 0, and system © can be written as 

. . /— y — asin(i) x + acos(t)\ 

where V(x,y,t) = (x + acos(t)) 2 + (y + asin(t)) 2 . This vector field corresponds to the 
following Hamiltonian 

H(x,y,t) = --hxV(x,y,t), 

that it is actually integrable. 
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Lemma 6.1. Let us consider e (energy), the symplectic variable conjugate to t, and define the 
autonomous Hamiltonian H\{x, y, t, e) = H(x, y, t) + e, then we have that 



-x 2 -y 2 



H 2 (x,y,t) = V(x,y,t)e~ 

is a first integral of Hi, in involution and functionally independent. As a consequence, if a = 
b 7^ the system is completely integrable. 

Proof. It is straightforward to see that the Poisson bracket with respect to the canonical form 
dx A dy + dt A de satisfies {Hi, H2} = 0. Moreover, 7i 2 does not depend on e, so that it is an 
independent first integral. □ 

Taking these results into account, one can completely understand the picture presented in 
the top-left plot of Fig. [3j Since the system is integrable, it is foliated by invariant tori, despite 
the two periodic orbits that are created by a resonance introduced when parameter A changes 
from A = to A ^ 0. Next, we characterize these two periodic orbits: 

Lemma 6.2. If A > 0, the system has two periodic orbits given by 

r ±(f) — ( x ±(t)yy±(f)) — ( a ± cost, a± sin t), 



where the coefficients a + and a_ are given by a± = ~ a± ^ 2+i . Moreover, the orbit r_ (t) is 
hyperbolic with characteristic exponents ±(ai — l) 1 / 2 , andr + (t) is elliptic with characteristic 
exponents ±i(l — a+) 1 ^ 2 - If A < the same result holds just switching the roles ofa + and a_. 

Proof. It is known that if the sets {H2 1 (c), c G M} are bounded differentiable submanifolds, 
their connected components carry quasi-periodic dynamics. Moreover, the critical points of 
H2 determine the periodic orbits of the system. Therefore, these periodic orbits are given by 
expressions: 2xV = and 2yV = which can also be written as 

x ((x + a cos t) 2 + (y + a sin t) 2 ) = x + a cos t, 
y ((x + a cost) 2 + (y + asint) 2 ) = y + asint. 

from which we obtain our two periodic orbits: x±(t) = a± cost, and y±(t) = a± sint, with 
a± = 1/ (a± + a). In addition, it is easy to check that a 2 _ > 1 and a 2 + < 1, respectively. 

Finally, the stability of these orbits can be obtained by considering the following associated 
variational equations, 

^ \ _ 2 ( sin(2t) -cos(2t) \( Wl \ 
w 2 J ± V -«>s(2t) -sin(2t) ) \ w 2 ) ' K ' 

Solutions for this equation can be easily obtained by using the complex variable z = w\ + 1W2, 
and solving z = —ie 2lt a 2 ± z. We have the following set of fundamental solutions 

Wl {t) = e^v 7 ^^ - a 2 ) cost T yj 'at - lsint^ , 
W2 {t) = e^V^^l - a 2 ) sint ± yj 'at - lcost^ , 
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for the hyperbolic case, and 



Wi(t) = cos(±tw 1 — a+)(±W 1 — a\ cost — (1 + a + ) sint) 



sin(±tW 1 - a%){±Jl - a% cost + (1 + a+) sint) 



iU2(t) — cos(±ty 1 — a+)(±y 1 — a+ sint + (1 + a+) cost) 

+ sin(±ty / ' 1 — a\){±\J 1 — a+ sint — (1 + a 2 + ) cost), 

for the elliptic one. Finally, the corresponding characteristic exponents can be obtained by a 
straightforward computation of the monodromy matrix. □ 

Remark 6.3. Notice that the chaotic sea observed in Fig. \3\ is associated to the intersection of 
the invariant manifolds of the hyperbolic periodic orbit that we have computed. 

Now, and in order to apply KAM theorem, we compute locally the frequency map of this 
unperturbed system around the vortex and the elliptic periodic orbit. To this end, we perform a 
symplectic change of coordinates in a neighborhood of these objects in order to obtain action- 
angle variables up to third order in the action. 

In general, let 7i(x, y, t, e) = H (x, y, t) +e be a Hamiltonian that is 2n -periodic with respect 
to t and has a first integral, F(x, y, t). Let us consider the generating function, S(x, t, /, E) = 
tE + S(x, t, /), determining a symplectic change of variables (x, y, t, e) t— > (J, 9, t, E) defined 
implicitly by 

- — e- E + — 9- — 
^ dx ' dt ' dl 

where 9 is also 2n -periodic. This transformation is introduced in such a way that the new 
Hamiltonian depends only on / 

Since a first integral of the system is known, we can define the corresponding action as 

I = F{x,y,t) = F(x,^,tj. (13) 

From Eq. (fT3l , we obtain locally the equation || = f(x,t,I), so that we have S(x,t,I) = 
J f(x, t, I)dx + g(t, I). Introducing this expression into (TT2l) . we obtain the following equation 
for g: 

^L = h(I)-H(x,f,t)- J^dx (14) 
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and can conclude that, since g must be 27r-periodic with respect to t, then h(I) has to satisfy 

h(I) = (H(x, /,*) + / %dx\ = (H(x, /, t)> , (15) 

where (■) denotes average with respect to t. Finally, we notice that since F is a first integral, we 
can define g so that 9 becomes 2n -periodic. 

Computations are simplified observing that the left hand side in Eq. (TBI) does not depend 
on x (we use the fact that F is a first integral), so we can set x = 0. According to this, we 
have to solve F(f, t) = I, where F(-, t) = F(0, •, t) and then we have to compute the average 
h(I) = (H(f, t)), where H(; t) = H(0, •, t). 

First, let us consider a neighborhood of the vortex for A > 0. To this end, we introduce 
the new variables x = —a cost + A x and y = —a sint + A y , so that the Hamiltonian H v = 
H V (A X , A y , t) and the first integral F v = F V (A X , A y , t) are 

H v = -^ln(A^ + Aj) - aA x cost - aA y sin t, (16) 



2 

F v = (A 2 X + A 2 ) exp(-a 2 + 2aA x cost + 2aA y sint - A 2 X - A 



Proposition 6.4. There exist a symplectic change of variables (A x , A y , t, e) i— > (I, 9, t, E), with 
9 G T, setting the vortex at I = 0, swc/z t/zat t/ze new Hamiltonian becomes 

h v {I) = -- In/ - - - — / + 3 (I). 

Proof. According to the above discussion, we have F v (f v ,t) = j2 e -a 2 +2a/„smt-/2 _ j j^ en ^ 
by introducing this expression in (fT6l) ones obtains 

- - 1 a 2 f 2 

H v (J v ,t) = --lnI---± 

Finally, we only have to compute the first terms in the expansion of f 2 obtaining 

f 2 = e a2 I- 2ael a2 sin t J 3 / 2 + 6a 2 e 2a2 sin 2 1 1 2 + . . . 

and use that (sint) = and (sin 2 1) = |. □ 

On the other hand, a neighborhood of the elliptic periodic orbit for A > can be studied 
by means of the variables x = a + cost + A x and y = a + sint + A r One thus obtains that the 
Hamiltonian and the first integral are given by 

H + = — \nV + + a + A x cost + a + A y sint, 
F + = V + exp(-a 2 + - 2a + A x cost - 2a+A y sint - A 2 - A 2 y ), 
where V + = (a + a + ) 2 (l + 2a + A x cost + 2a + A y smt + a^A 2 + a+A 2 ). 
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Proposition 6.5. There exist a symplectic change of variables (A x , A y , t, e) i— > (I, 8, t, E), with 
9 e T, setting the periodic orbit at I = (a + + a) 2 e~ a , such that the new Hamiltonian becomes 

h + (I) = - ln(a + a + f + J - 1 -±^kj* + Q 3 (I), 

where we have introduced the notation 

e a2 I 



J=l 



(a+ + a) 2 
and also 

a 2 (41a 8 - 88a 6 + 119a 4 - 54a 2 + U 



ni = r. — r , n 



y/1 - a 4 ' 3QVT^a I (a 8 + 1 - 2a 4 )(l + a 2 ) 

Proof. As before, we consider a solution /+(/, t) for the equation F + (f + , t) = I. For conve- 
nience, we introduce the notation / = (a + + a) 2 e" a (1 — J) in order to set the periodic orbit at 
J = 0. Then, it turns out that the expression 

(1 - a\ cos(2t))/ 2 + [a + {\ + a 2 ) sin t - j|a 3 sin 3 tj / 3 + 4 (/+) = J, 

approximates the previous equation for /+ and that the following expansion in terms of J 

fl = ai {t)J + a 3/2 (t)J 3/2 + a 2 (t)J 2 + ... 



holds, where 



1 

OL\(t) 



«3/2(» 



1 — ai cos 2t 

+ 



a+(l + a 2 ) suit + fa 3 sin 3 t 



[1-al cos (2t)) 5 / 2 



_3(a+(l + a+) 2 sint- §a 3 sin 3 t) 2 
Q2{t) = 2 (l-a 2 + c os2t) 4 ' 

Hence, we have to compute the average of 

H+(f+,t) = - In(o+ + a) 2 - 1 ln(l - J) - 5 

that follows from the fact that (cui) = III, (a 3 / 2 ) = and (a 2 ) = n 2 . These averages are 
computed easily by using the method of residues. □ 



16 



Bohmian trajectories from a Dynamical Systems point of view 





Integrable 


Hamiltonian 


Time 
Reversible 


Stroboscopic 
sections 


A = 0,B ^0,C 7^0 


yes 


no 


yes 


Ellipses 










around origin 


A^Q,B = C 


yes 


yes 


yes 
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no 


no 


yes 


Rest of panels 










in Fig. [3] 



Table 2: Dynamical characteristics of the quantum trajectories generated from the different possibilities in the 
canonical model ([H) for the pilot wave function ©. 

7 Conclusion 

In this paper we present an scheme to study in a systematic way the intrinsic stochasticity and 
general complexity of the quantum trajectories that are the basis of quantum mechanics in the 
formalism developed by Bohm in the 1950's. In our opinion this approach, which based on 
the ideas and results of the dynamical systems theory, can seriously contribute to establish firm 
grounds that foster the importance of the conclusions of future studies relying on such trajec- 
tories, thus avoiding errors and ambiguities that has happened in the past. As an illustration 
we have considered the simplest, non-trivial combination of eigenstates of the two dimensional 
isotropic harmonic oscillator. 

The corresponding velocity field is put in a so-called canonical form, and the characteristics 
of the corresponding quantum trajectories studied in detail. It is proved that only one vortex 
and two periodic orbits, one elliptic and the other hyperbolic, organize the full dynamics of 
the system. In it, there exist invariant tori associated to the vortex and the elliptic periodic 
orbit. Moreover, there is a chaotic sea associated to the hyperbolic periodic orbit. The KAM 
theory has been applied to this scenario by resorting to a suitable time-reversible symmetry, 
that is directly observed in the canonical form for the velocity field determining the quantum 
trajectories of the system. It should be remarked that the results reported here concerning the 
hyperbolic periodic orbit constitute a generalization of those previously reported in |[T5l . in 
the sense that here a more concise and constructive approach to the associated dynamics, is 
presented. We summarize the dynamical characteristics of the different possibilities arising 
from the canonical velocity field ([8]) in Table [21 that represents a true road-map to navigate 
across the dynamical system, i.e. quantum trajectories, that are defined based on the pilot effect 
of the wave function ©. Also, note that the generic model, i.e. when E, F or G do not 
vanish, does not satisfy any of the properties considered in the table. 

Finally, the method presented here is, in principle, generalizable to other more complicated 
situations in which more vortex and effective dimensions exist. Some methods have been de- 
scribed in the literature that can be applied to these situations |[30ll . This will be the subject of 
future work. 
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